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Nonequilibrium transport properties are determined ex- 
actly for an adiabatically contacted single-channel quantum 
wire containing one impurity. Employing the Luttinger liquid 
model with interaction parameter g, for very strong interac- 
tions, g < 0.2, and sufficiently low temperatures, we find an 
S-shaped current-voltage relation. The unstable branch with 
negative differential conductance gives rise to current oscilla- 
tions and hysteretic effects. These nonperturbative and non- 
linear features appear only out of equilibrium. 

PACS: 72.10.-d, 73.40. Gk 



Transport in ID conductors is one of the focal points 
of condensed matter physics. At low energy scales, such 
materials have been predicted long ago to behave as Lut- 
tinger liquids (LL) instead of Fermi liquids EL Over the 
past few years, several possible experimental realizations 
of LL behavior have been reported. In particular, narrow 
quantum wires (QW) in semiconductor heterostructures 
can be operated in the single-channel limit [|]||]. Other 
realizations include quasi-lD materials such as long chain 
molecules [0 , carbon nanotubes or edge states in frac- 
tional quantum Hall (FQH) bars iBj. Since the latter are 
in fact chiral LLs, where right- and left-moving branches 
are spatially separated, FQH edge state transport j7|,|| is 
distinct from the case of a quantum wire. In this Letter, 
we emphasize the important and indeed surprising differ- 
ences arising for standard (achiral) LL systems charac- 
terized by the interaction parameter g < 1. We focus on 
the archetype problem of a spinless single-channel QW 
containing backscattering (BS) by one impurity 0. 

Our main results are as follows. The current I(U, T, g) 
under an applied voltage bias U at temperature T obeys 
scaling [I depends only on U and T measured in terms of 
the impurity scale Tg] and a duality relation connecting 
the strong and weak BS limits under the simultaneous 
exchange g — > Both the scaling function and the 

duality relation are different from the FQH case and are 
determined exactly. For very strong interactions, g < 0.2, 
and low temperatures, the I — U characteristics is multi- 
valued, containing an unstable branch of negative differ- 
ential conductance (NDC). Once the QW is embedded 
in a load circuit, this S-shaped I(U) relation can lead to 
hysteresis, current switching and self-sustained current 
oscillations ||. Such effects could be observed in a QW 
of very low electron density. 



Coupling of QW to voltage bias. — We focus on a 
QW adiabatically connected to ideal time-independent 
voltage reservoirs held at electro-chemical potentials /L«i,2 
The ID conductor extending from -L/2 < 
x < L/2 is considered to be a LL, with an impu- 
rity of BS strength A sitting at x = 0. This defines 
T B = CgXiX/tOcy/^-a^ where % = k B = 1, oj c is the 
electronic bandwidth, and c g is a numerical prefactor of 
order unity (its precise value is of no interest here and 
given in Ref. jt3|). The Coulomb interactions take the 
form Hi = (e/2) J dxp(x)ip(x), where p(x) is the electron 
density and the Poisson equation is replaced by 

eip(x) = uop(x) with g = (1 + uq/ttvf)^ 1 ^ 2 , (1) 

with the Fermi velocity vp. Here a screening backgate or 
other metallic surroundings cause short-ranged interac- 
tions within the QW. Since we are dealing with a strongly 
correlated non-Fermi liquid system, the well-known Lan- 
dauer approach |lC| ] does not apply. In the past, external 
voltage sources were often modeled by attaching g = 1 
LLs to the ends of the QW |Q , but computations become 
exceedingly difficult for A > 0. Therefore we employ the 
radiative boundary conditions of Ref. pd| ] which repre- 
sent the natural extension of Landauer's original ideas 
|l0f to a strongly correlated ID metal. Using different ar- 
guments, these boundary conditions have been confirmed 
and generalized to a.c. transport |l^ , [l7|] . 

Let us briefly summarize the main ideas fll]]l2"|| . Sup- 
pose one injects the "bare" densities p R and p\ at posi- 
tions x close to the end of the QW. The average density 
p = Pr + Pl [we omit the expectation values for brevity] 
is then self-consistently determined by 

p R (x) + p L {x) = p% +p° L - eLp(x)/irv F , (2) 

since the band bottom shifts by eip(x). With Eq. (0) we 
then get the local relation pn + pl = 9 2 (Pl + P%), in 
agreement with the compressibility of a LL, k — g 2 /ttvf- 
Since screening affects only the total charge, see Eq. (Q), 
we also have pr — pl = p% — p\- Solving these two 
relations for p R j L and using p R (—L/2) = p\j2-nvF and 
p L {L/2) = /J2/27TUF, we obtain the boundary conditions 
ofRefs. OMR!, 
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where we put p\ = —pi = J7/2 with the applied voltage 
U and e = 1. We assume full translational invariancc 
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such that the sound velocity v = vp/g (below we set 
v = 1). The Sommerfeld-like boundary conditions (||) 
are imposed at the left/right end of the conductor at 
long times t where the stationary nonequilibrium state 
has been reached. Below we assume that the relevant 
energy scale (fcsT or elf) exceeds v/L, and henceforth 
take L — > oo. 

Exact solution and duality. — Like for the tunneling 
problem in the FQH effect M, it is possible to compute 
exactly the current out of equilibrium. The basic idea is 
first to fold the problem onto the boundary sine-Gordon 
model, and then use integrability of the latter |l3[ . The 
final equations governing the physics are quite simple. 
For clarity, we restrict ourselves to the case g = l/p 
with p integer. The basic quantities are then pseudo- 
energies £j(0) for rapidity 9, obeying a set of thermody- 
namic Bcthe ansatz (TBA) integral equations, 
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where s(9) — (p — l)/cosh[(p — 1)9} and Njk is the in- 
cidence matrix of the following TBA diagram, on which 
the labels j, k run: 
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Here we have defined parameters for the breathers m. 



2 sin ■ 



with j 



1, . . . ,p — 2, and for the kink and 
. The latter are the fundamental 



_ ■ 2(p-l. 

antikink, m± = 1 
charged particles in the integrable description. All the 
physical particles arc coupled by non trivial S matrix 
elements. However, manipulation of the TBA equations 
gives rise to the simpler structure of coupling for the pseu- 
doenergies represneted on the foregoing diagram. The 
chemical potentials are [ij = for the p— 2 breathers, and 
/i ± = for kink and antikink. Here W is determined 

self-consistently, see below. Having found the e's, the 
densities of kinks and antikinks are given by <r± — nf± , 
where the pseudo-energies e± are equal, 27m = de±/d9, 
and the filling fractions read 



f±(9) = 1/{1 + exp[(e(0) T W/2)/T}} 



(5) 



With these definitions, the final expression for the current 
reads 



(0) 



1= / \T++\'(a+-a-)d0, 



where, with the impurity scale Tb oh exp 9b defined 
above, the tunneling probability is 



|T ++ | 2 = (l + cxp[-2( g - 1 -l) (9 -9 B )] 



(7) 



The current is implicitly a function of W, which is self- 
consistently determined through 

/ V + _i 2 )(a + -o^ = |;, (8) 

where |T+_| 2 = 1 - |T ++ | 2 . 

^From these equations, it is now straightforward to 
deduce the following identity giving the parameter W in 
terms of the physical voltage and current, 



U = 2tv[1--]I + W 

9, 



In the sequel, we will also use the quantity 
V = U - 2-kI = W - 2nl/g , 



(9) 



(10) 



whose physical meaning is the four-terminal voltage 
across the impurity H], i.e., the voltage difference mea- 
sured by weakly coupled reservoirs on either side of the 
impurity. We then find that for non- vanishing Tg, the 
current interpolates between 1 = and (going back to 
physical units) I = (e 2 /h)U as the applied voltage U is 
increased. Similarly, the linear conductance interpolates 
between the perfectly quantized high-temperature value 
G = e 2 /h also found in a clean QW jl4j], and the vanish- 
ing low-temperature (T <C Tb) conductance predicted in 
Ref. 0. For temperatures above ui c or vp/r, where r is 
the interaction range, the conductance saturates before 
reaching e 2 /h in a real QW. 

The above equations can be solved in closed form at 
vanishing temperature. The solution for arbitrary g is 
expressed in terms of two different series expansions, de- 
pending on whether the impurity BS is weak or strong. 
In the latter case, we find 



I = G(p)(e A Mj2(-l) n+1 



IH T(np) 



2T(n)T(n(p- 1) + 3/2) 



o A+A-0i3\ 2 ™(p-l) 



x (e— , (11) 

while the boundary condition (^|) reads 

oo 

U = 2G(p)e A - (p - l)G(p)e A ^(-1)" +1 

71 = 1 



X ■ 



r(n)r (n(p - 1) + 3/2) 



(12) 



Here we have introduced the notations G{p) = 
p^T{p/2(p - 1)) / r(l/2(p - 1)) and A = [ln(p - 1) - 
p\n(p)/(p — l)]/2. The parameter A follows by solving 
Eq. ( |l2| ) for a given voltage U . Inserting it into Eq. ( pd| ) 
gives the I — U relation and W = 2G(p)e A . Equations 
in the weak BS limit follow from the duality relation de- 
scribed below. 
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The TBA equations are easily solved for any T at g 



1/2, with the result (12 



(13) 



where ip is the digamma function and Tb = ir\ 2 /lu c . The 
linear conductance is then 



e 2 1 
G(T) = 6 - 



h i + c 0/(i+c) 



c = T b /2ttT , (14) 



with the trigamma function ij)'(x). For high tempera- 
tures, this approaches e 2 /h, while it vanishes as G ~ T 2 
at low temperatures. Notably, while this is the same 
power law as in the FQH effect fij, the prefactor is now 
different. To get the current for arbitrary values of g,T 
and f, one has to resort to a straightforward numerical 
solution of the TBA equations ( |ll| ) and ([12]). The linear 
conductance can be given in closed form, 



G = 



f de\T ++ \ 2 df±/de 



f d6(\T ++ \ 2 + \T + 4 2 /g)df ± /d6 



(15) 



where the filling fractions (||) have to be evaluated at 
W = 0. 

A remarkable nonperturbative consequence of the 
TBA equations is the existence of a duality relation for 
the current valid at any temperature, 



I(X,U,g) = ^--l(x d ,U,- 
2tt V 9 



(16) 



where A^ oc A -1 / 9 jlj|. This duality is similar but differ- 
ent from the one found in the FQH case Although, 
like in the FQH case, quasiparticles of charge q = ge 
tunnel in the weak BS limit, and electrons with q = e 
in the strong BS limit, the coupling of the QW to the 
reservoirs modifies the current and gives rise to the new 
duality @. 

Bistability regime. — The curves for the linear con- 
ductance in a QW are generally similar to those for the 
FQH effect, with the main difference that the high tem- 
perature value is now e 2 jh independent of g. Much more 
interesting and unexpected physics occurs in the nonlin- 
ear out-of-equilibrium regime for very strong interactions 
and low temperatures. The most direct approach to see 
this is the quasiclassical limit, g <C 1, at T = 0, which we 
consider first. In the bosonized theory, after integrating 
out the standard boson phase fields away from x = but 
taking into account the boundary conditions (||) , one 
is left with the equation of motion 



d^/dt + q-kTb sin $ = gW , 



(17) 



where Tb — 2A and the current operator is <f>/2ir. This 
equation can also be obtained from the explicit solution 



of the TBA in that limit. For W > nT B , the current 
grows proportional to A = [W 2 — (ttTb) 2 ] 1 ^ 2 . This is 
readily seen from the solution of Eq. (|l7j), which reads in 
terms of y — tan [g At/2], 

tan(*/2) = Wy(t)/[A + irT B y(t)} . 

Hence the average over one time period yields 



(/) = 7^— A Q(W — ttTb) 



(18) 



where O is the Heaviside function. The I(W) curve is 
single-valued, but by inserting the definitions (Q) and 
(|Io|), and eliminating 7, we see that U — F(W) with 



F(W) =gW+(l- g)V{W) 



(19) 



can have several solutions W . Therefore the I(U) rela- 
tion is multivalued throughout the regime U\ < U < U2, 
where U\ = Wg(2 — g) ttTb and U% = ttTb- The three 
solutions for the current in this regime are 1 = and 

U 



I±(U) = 2n{2 _ g) (l - 9 ± V~l - g (2 - g)(*T B /Uy 



(20) 



For U < Ui, we only get 1 = 0, and for U > U 2 , I+(U) is 
the only allowed solution. Clearly, on the branch 7_(Z7), 
we have negative differential conductance (NDC), and 
therefore this branch is unstable. Such S-shaped current- 
voltage relations are familiar, e.g., in nonlinear semicon- 
ductor physics H and in charge-density wave transport 
|l9| . By putting the QW into a properly designed load 
circuit, self-sustained current oscillations can be gener- 
ated Q. Furthermore, putting a resistance R in series 
and applying the voltage V a to the whole circuit, one has 



V a = RI(U) + U 



(21) 



which can easily be solved for I(V a ). For R/(e 2 /h) > 
g^ 1 — 1, we get a single- valued I(V a ) curve, i.e., the NDC 
branch has been stabilized. For smaller R, also the I{V a ) 
curve will be multi- valued. In practice, one then gets 
hysteresis, and the current jumps between two branches 
(bistability) . 

How stable is this behavior once thermal and quantum 
fluctuations are taken into account? Repeating the above 
g -C 1 calculation for finite temperature with methods 
from Ref . pQl , one finds 



V(W) = ttT b Im — 

1 -iW/2irT\ 1 B/ £± ) 



(22) 



with the modified Bessel function I v (z) with complex or- 
der v. The high-temperature expansion (T ^> Tb) of 
Eq. ( p2| ) predicts bistability to occur for all T < T c (g) 
with 
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T c (g<Zl)=T By /(l-g)/16g. (23) 

Therefore as g — > 0, the bistability occurs at all temper- 
atures. Bistability is not destroyed by quantum fluctua- 
tions either, as can be checked using the exact solution. 
At T = 0, one finds that the NDC disappears for g w 0.2, 
see Figure |l[ This value gets smoothly lowered as tem- 
perature increases. We stress that the bistability is a 
true nonperturbative nonequilibrium phenomenon that 
has no parallel in the linear conductance pj ]. 

Finally we turn to rather elementary physical reason- 
ing explaining the origin of the predicted bistable behav- 
ior independent of the details of our boundary condition. 
In the presence of impurity BS, a portion V of the volt- 
age drops at the impurity site. Then the average current 
(I) = (U — V)/2n, where l/2n [— e 2 /h in physical units] 
is the d.c. conductance of the perfect wire. However, a 
fluctuation SI of the current must obey 

SI = -{g/2n)SV , (24) 

with the associated fluctuation SV of the four-terminal 
voltage and the a.c. conductance g/2ir 0. Equation ( |2~i| ) 
reflects the fact that the impurity BS potential Acos$ 
for the x = boson field <& is due to tunneling of frac- 
tional quasiparticles with charge ge p2[ . In terms of 
the boson field, the current is (/) + SI = $/27r, and 
the voltage fluctuations read 5V = 2ir\ sin <I> — V. The 
first term is basically the derivative of the pinning po- 
tential Acos$ p2| . With these relations, it is straight- 
forward to verify the equation of motion ([l7]). The lat- 
ter describes an overdamped particle moving in the tilted 
washboard potential — g(2nX cos & + where the bias 
W = U + 2ir(g~ 1 — l)(I) depends explicitly on the cur- 
rent. This feedback mechanism together with the nonlin- 
ear pinning potential is responsible for bistability. The 
above arguments also suggest that our assumption of adi- 
abatic coupling to the reservoirs is not essential for bista- 
bility to occur |^3[ . 

To conclude, nonequilibrium transport through a spin- 
less single-channel quantum wire containing one impurity 
has been studied. Using integrability techniques, the ex- 
act solution of this interacting transport problem has 
been given for adiabatically connected reservoirs. We 
have discovered bistability phenomena in the current for 
strong interactions that should be observable in state-of- 
the-art experiments and constitute a hallmark of strongly 
interacting quantum wires. 
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FIG. 1. Current-voltage relation at T = for various 
g = 1/p. The current is given in units of (e 2 //i)Ts, and the 
voltage in units of Tb ■ 
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